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ABSTRACT 

This paper explores fluctuations and noise in various facets of cancer development. The three areas of particular 
focus are the stochastic progression of cells to cancer, fluctuations of the tumor size during treatment, and noise 
in cancer cell signalling. We explore the stochastic dynamics of tumor growth and response to treatment using 
a Markov model, and fluctutions in tumor size in response to treatment using partial differential equations. 
We also explore noise within gene networks in cancer cells, and noise in inter-cell signalling. 
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1. INTRODUCTION 

Fluctuations and noise are present in a wide variety of physical systems. This paper explores the effects 
of noise in tumor growth and response to treatment, and in general "noise" in gene networks and in inter-cell 
signaling. This improves our understanding of the nature of cancer development and growth, through its 
examination of the robustness and sensitivity of the gene pathways that are important in the progression of 
normal cells to cancerous cells. 

2. METHODS 

For both the tumor response and Wnt signaling analysis, the ODEs are solved using the Runge-Kutta method 
of order 5 with variable step size. We note that although an Euler method has been shown to work better than 
the Runge-Kutta method, it has been shown that using the Runge-Kutta method of order 4 with step sizes 
much greater than those used here (of a maximum h — 0.1) works well for both the linear additive noise used 
in the Wnt signaling modeling and in particular for the multiplicative noise used in the tumor response.^ 

We explore three models useful for studying cancer. Two models consider the effect of tumor treatment, 
one a Markov model and the other an ODE model. In the Markov, population-based model, we consider 
stochastic effects on tumor growth when treatment is applied. Observing the fluctuations in the response leads 
us to consider adding such fluctuations to an existing ODE model of tumor treatment,^ which considers the 
cells in regions of 3D space. We also explore a model of inter-cell signalling and the effect of adding noise to 
it, as inter-cell signaling plays an important role in the development of cancer.® 

Send correspondence to Derek Abbott 
E-mail: dabbott@eleceng.adelaide.edu.au, Telephone: -1-61 8 8303 5748 



2.1. Tumor treatment using Markov model 

We consider in our models the following cell populations: normal epithelial eells denoted by N eells with 
mutations that lead to genetic instability (G), cells with mutations which increase their replication rate (R), 
cells with mutations which allow them to avoid death (£>), cells which have acquired the ability to induce 
angiogcnesis (^4). and cells with two or more mutations. Cell populations that have acqiiired two or three 
mutations are denoted by listing the mutations together. We label a cell which has acquired all four mutations 
a tumor cell (T). Finally, a tumor cell which has acquired the capability to invade and metastasize is a cancer 
cell {M). 

The spontaneous mutation rate in human cells has been estimated to be in the range of 10~^ to 10~^ 
mutations /gene /cell division.^ We assume a mutation rate of ki = 10~^/gene/ccll division. As genetic 
instability favors the accumulation of mutations,* we assume that the mutation rate after a genetic instability 
mutation increases to k2 = 10~^ mutations/gene/cell division.^ Successful invasion and metastasis depend 
upon acquisition of the other capabilities, as well as several new capabilities.^" To simplify the model, we do 
not address the multistep progression of a tumor cell to a metastatic cell, and instead consider this complex 
process as one step. It is estimated that only 1 in 10^ cancer cells will invade and successfully metastasize to 
a new site,^^ and so we use a rate of ^3 = 10^^ for the transition from a tumor cell to a metastatic cell. In 
addition, our model assumes that a tumor cannot grow past « 10'' cells without angiogenesis supplying blood 
to the tumor, and thus we cap the size of the tumor at Ci — 10® cells until the tumor acquires a mutation 
in an A gene. Death of the patient occurs when the tumor reaches C2 = 10^^ cells. The Markov model also 
takes into account the fact that there are approximately 100 genes involved in each category of mutations.^ 

The observed volume doubling time of tumors in clinical practice is on the order of 100 days."'^^ We assume 
that the relative contribution to increased net proliferation for mutations in the D and R categories is 0.8 and 
0.2, respectively. This ratio places more importance on apoptosis mutations than the 0.7 : 0.3 inferred from 
a paper by Tomlinson and Bodmer^''^ due to the crucial role of apoptosis in net cell proliferation.^*' There 
is no net proliferation associated with cells which do not have mutations in D or R. The net doubling times 
associated with D and R in our models are 500 and 200 days, respectively. Cells with a mutation in R and D 
have a net doubling time of 100 days. Cells without a mutation in an R gene divide every 6 = 5 days; cells 
with a mutation in an R gene divide every bn = 4.96 days. The cell lifetime used for cells without a mutation 
in a f gene is d = 5 days; the lifetime used for cells with a mutation in a D gene is do = 5.08 days. The birth 
and death rates arc equal for normal cells, thus tlica'c is no net proliferation unless a mutation in D or R is 
acquired. The cell division rates are conservative estimates based on Rew and Wilson.^'* 

Our model makes several simplifications. In the case where some cells might not have a mutation in A but 
neighboring cells do, we still limit the population of cells without mutations in A even though they may have 
an adequate blood supply due to signals from neighboring cells. In addition, we do not take into account the 
fact that mutations in some genes may fall in two or more categories. For example, a mutation in p53 may 
fall into both the D and G categories, as p53 is involved in both apoptotic and DNA repair pathways. 

We consider a discrete time Markov model, using the states and transition rates as shown in Figure 1 . We 
use state-dependent time steps of At = 1/6, where 6 is the cell division time, as we consider the mutations 
to become fixed at each cell division. For example, if a cell is in the R state, each bn days, the cells gain 
mutations with probability 

P=Ps{t)k,g, (1) 

where Ps{t) is the population of cells in state s at time t, kg is the rate of transition out of state s, and g is 
the number of genes involved in each state, which we keep constant at 100. We begin with a population of 
normal cells (except when we are modeling inherited mutations) . We then allow the cells to mutate through a 
particular path of states, dividing at a rate which depends on the current state. The net proliferation depends 
on the current state as well. We treat the total population of cells as a homogeneous population, and run 
the simulation for a large number of trials (10,000). We then average in order to simulate a heterogeneous 
population. 




Figure 1. State diagram of tlie model. Normal cells (A'') can acquire mutations which lead to genetic instability, (G); 
mutations which increase the proliferation rate, {R); mutations which give the cell the capability to avoid death, {D); 
or mutations which give the cell the capability to induce angiogenesis, [A). These mutations are acquired at rate ki. 

After a mutation in G, the mutation rate increases to k-z. When a cell has acquired a mutation in each category, it 
becomes a tumor cell (T). Finally, tumor cells become metastatic cells (M) at rate fcs. 

2.2. Tumor treatment using an ODE model 

Wc use a partial differential equation l)ased model for tumor size based on work by Jackson and Byrne."'' Tire 
model treats the tumor as a spherical mass with rapidly dividing, highly drug-susceptible cells, surrounded by 
normal cells that have a lower response to the drug treatment. The parameters used are defined in Table 1, 
and the symbols in Table 2, as used in Equations 2-4. 

Table 1. List of the parameters used in the model and the values used from the literature.^' ^® The parameter F 
is the rate coefficient of blood-tissue transfer, or the rate at which the drug can pass from the blood into cells. The 
parameter is the association rate, or binding of drug to target receptors in the tumor. The rates of drug loss due to 
decay, molecular instability, or cellular uptake and metabolism are combined in the parameter A. Transfer coefficients 
ki2 and k2i are used to describe transfer of the drug out of and into the group of rapidly dividing, highly responsive 
cells, respectively. Parameter ke represents elimination of the drug from blood plasma. 



parameter 


value 


r 


16 day-^ 


/3 


560 M-May-^ 


A 


1.9 day-^ 


h 


2.4 X 10-=^ s-i 




5.8 X 10"'' s-^ 




7.1 X 10-'^ s-^ 



Table 2. Variables used in the model. 



symbol 


description 


du 


jjicscrilKHl drug coiicculrMtion in the tumor vasculatmc (M) 




the drug concentration in the surrounding tissue (M) 


R 


the radius of the tumor (mm) 



-TT = —ki2dB + k^idN — keds, 



(2) 



^ = ki2dB - k2idN, (3) 
dt 

dR / /3r, ..\R (3(, r, , , \ ^i?coshCi? - sinh^i? 



where fee, A:i2, and k2i are given in Table 1, and ^ = VT + A. The initial condition for d^, C, is given by 

C = ^Jkl^ + 2ki2ke + 2/ci2fc2i +kl + fc^i/fcia, (5) 

where fci2, ^21, and /ce are as defined in Table f. The initial conditions used are i?(0) = 0.4 cm"^ (in 
line with Kerr et al}^), dB{0) = 0.1, and djv = C We consider fluctuations of the tumor size, R, in the 
response to the drug concentration in the tumor, ds- We multiply an exponential noise term a, from an 
exponential distribution P{x) = Xe^^^ , by ds to show fluctuations in tumor response in short intervals, and 
decreased tumor responsiveness to treatment over time. Tumor responsiveness may vary with fluctuations in 
cell signaling, which we explore in the following subsection. 

2.3. Inter-cell signalling 

The Wnt signaling pathway contains many oncogenes. Although the evidence for mutations in and ovcrex- 
pression or midcrcxpression of Wnt in development of cancer is not entirely conclusive, there is growing 
evidence that the downstream targets such as APC and /3-catenin play an important role in the development 
of many cancers.^' ^"'^""^ 

Here we consider the reaction scheme for Wnt signaling shown in Figure 2. In this reaction scheme, Wnt, 
through binding to its receptor Frizzled, activates Dsh. This activated Dsh protein then binds to axin, and 
through bound GBP, inhibits the phosphorylation of /3-catenin. Axin plays a key role in coordinating the 
assembly of large protein complexes including one that forms from glycogen synthase kinase 3/3 (GSK3/3), 
adenomatous polypopsis coli (APC) and the the negateivc regulators Dsh and GSK3/3- binding protein (GBP). 
The absence of Wnt signals for the phosphorylation of /3-catenin, marking it for subsequent degradation in the 
degradation cycle shown. Since /3-catenin plays a role in signaling events in the cell-cycle, the Wnt pathway 
plays a role in tumor formation, where the ceU-cycle becomes disrupted. 

The Wnt signaling pathway can be modeled using a set of ordinary differential equations. We use the 
same set of equations in modeling /3-catenin degredation, however we expand the range of functions used to 
model the transient stimulation of Wnt in vitro, and model various types of noise on these inputs. 

Although the degradation of Wnt is assumed to be an exponentially decaying function in Lee et al.,'^'^ this 
assumes that: 

• the degradation process is random, 

• independent of time, 

• has a constant probability per unit time for each Wnt molecule. 

It is not obvious biologically that the process is entirely random, nor that the probability per unit time is 

entirely constant. While clearly the end result of a decrease in Wnt signalling is an exponential-like decay in 
/3-catenin, there are many simpler expressions which also give such a decay in /3-catenin. 

The various functions for the Wnt degradation, W (t), we use are given in Equations 6 through 8. First, 
we consider the simplest function, that the Wnt is degraded linearly: 

W{t) = h^'^' (6) 
[max{l- X{t-to),0} + n{t), t > to, 




Figure 2. This figure (from Lee et al.^^) sliows tlie /3-degradation cycle, and tfie protein interactions are sliown with 
arrows. 



where A = 1/tw for tw the characteristic lifetime of 20 minutes for Wnt, to is the time at which the signahng 
occurs, and n{t) is the noise term described later. This linear model changes the assumption of constant 
probability per unit time for each molecule; instead it describes a constant rate of breakdown independent 
of the amount of Wnt. In other words, this represents saturation of Wnt receptors. We also tried the same 
exponential decay equation as used by Lee et al, however with a noise term, n{t), added: 

W(t)^f^' (7) 
|e-^(*-*°)+n(t) t>to, ^' 

We used a combination of Equations 6 and 7, to reflect saturation of Wnt receptors, as in Equation 6, to the 
point where it can be degraded in an independent fashion as in Equation 7: 

{n{t), t<to, 
1- X{t~tQ)+n{t), to<t<l/\ + Wc/X + to, (8) 

^^g- A(t-i/A+wjA) +n{t), t>l/X + Wc/X + to/X, 

Where < Wc < 1 is the level at which Wnt can be degraded in an independent fashion. The types of noise, 
n(t) we consider are: 

• uniform noise in the range (0,0.1), which represents a small, approximately constant background level 
of Wnt over time, 

• Poisson noise, representing a a level of Wnt which occasionally fluctuates to high levels through time, 
and 

• exponential noise, similar to the Poisson noise but with no structure to the time between fluctuations. 



3. RESULTS 
3.1. Tumor treatment using a Markov model 

Figure 3 shows the effect of applying surgery followed by intermittent therapy to tumor masses of various sizes. 
The stochastic nature of the Markov model is visible, although use of a large number of simulations (10, 000) 
reduces this somewhat. 

We show a noticeable drop in cell population in Figure 3(a) as therapy effectiveness, the percent of cells 
being killed on replication, increases from 70% to 90%. The minimum effectiveness to significantly reduce the 



Table 3. Here we show the level of /3-catenin (in nM) remaining 180 minutes after removal of steady-state Wnt signal, 
for all three types of decaying input signals and with no noise and all throe types of noise added. 





Exponential input 


Linear input 


Mixed input 


No noise 


9.48 


8.07 


9.15 


Exponential noise 


11.36 


10.79 


11.21 


Uniform noise 


11.05 


10.21 


10.84 


Poisson noise 


9.89 


8.43 


9.53 



population of metastatic cells is reasonably realistic here as the Markov model captures the discrete time nature 

of the intermittent treatment. The transition from very little effect to a significant effect on the population of 
metastatic cells is fairly gradual between 70% and 90%. The effect of this increase in treatment effectiveness 
on cell population is shown in Figures 3(b) and 3(c), where we plot the number of M cells over time with 
varying tumor sizes at the start of treatment. 

From Figures 3(b) and 3(c), it is clear that the treatment effectiveness is more important than the tumor 
size at which therapy is applied. With a therapy effectiveness of 70%, there are w 10^ M cells at 100 years on 
average, and with a therapy effectiveness of 90% there are no M cells at 100 years on average. 

The level of fluctuations in tumor size increase over time, highlighting the variation in tumor response 
to treatment as some cell populations within the heterogeneous tumor mass gain resistance to the treatment 
drugs used. This is one reason why we consider such variations in the tumor response ODE model. 

3.2. Tumor treatment using a DE model 

Normally, application of the treatment holds the growth of the tumor in check, as shown in Figure 4(a). The 
tumor expands to a size at which a large amount of the drug can enter, then starts to shrink reasonably slowly. 
Multiplying the tumor response to treatment term, dsit) by exponential noise, with probability distribution 
P{x) = Xe~^^, for a wide variation in A (0.5 through 1.5), produces uncontrolled tumor growth within a very 
short period of time, as shown in Figure 4(b). Our results highlight the importance of considering fluctuations 
in the tumor response. The difference they make to the model is stark: we have gone from a very controlled, 
good response to treatment to a situation where the tumor grows uncontrollably. Such fluctuations may be 
caused by variations in cell division rate (as shown by the Markov model) which affects the response to the 
drug, or other cellular responses such as inter- and intra-cell signaling. 

3.3. Inter-cell signalling 

The results of using the linear, c;xponential, and mixed rates of Wnt decay as per Equations 6-8, and with no 
noise added, are shown in Figure 5. All the various forms of noise added to the system increased the level of 
/3-catenin at 180 minutes, with exponential noise increasing it the most. The values are shown in Table 3. For 
uniform and exponential (correlated) noise there is a higher overall background level of Wnt, which translates 
into a higher level of /3-catenin in the system, when compared with Poisson (uncorrelated) noise. We then 
added in the set of noise terms, n{t) into Equations 6-8. The result of adding in noise is shown in Figure 6. 



0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 

Therapy effectiveness 



(a) 




120 



Time (years) Time (years) 

(b) (c) 

Figure 3. Simulations of therapy using the Markov model with varying parameters of effectiveness and cell population, 
(a) Number of metastatic cells at 100 years after applying surgery at 10^ tumor cells followed by chemotherapy with 
varying degrees of effectiveness, (b) Effect of therapy that is 70% effective, applied at different tumor sizes: 10*, 10^, 
and 10^" . (c) Effect of therapy that is 90% effective, applied at different tumor sizes. 
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(a) Growth of tumor with treatment, no fluc- 
tuations in response. 



(b) Growth of tumor with treatment, fluctua- 
tions in response. 



Figure 4. ODE model of tumor growth and response. This figure demonstrates the controlled tumor growth for the 
case where treatment is applied with no fluctuations in tumor response (a) and uncontrolled tumor growth for treatment 
with minor fluctuations (b). In (a), the tumor initially grows in size up to the point where the surface area is large 
enough that the drug can enter the tumor in high enough concentrations to start having an effect. Note the difference 
in scale of the y-axis between the two graphs. 




— exponential input 
poisson input 

- - uniform input 



(a) This shows the linear, mixed, and exponen- 
tial degradation of the Wnt input signal over 
time, with no noise added. 



(b) Here wc show the response of the system 
for varying inputs. We plot the concentration 
of /3-catcnin over time for the exponential, lin- 
ear, and mixed linear and exponential decrease 
in the level of Wnt. 



Figure 5. Input and output of the Wnt signaling system in the absence of noise. Note that all inputs give exponential 
decays of /3-catenin, however a linear decrease in Wnt causes a larger decrease in /3-catenin. This is due to the fact that 

for the linear decrease, the Wnt receptors arc being saturated and thus signal for degradation of /3-catenin at a higher 
rate. The mixed model appears somewhere between the two, as the Wnt receptors are only initially saturated. 




Figure 6. These figures show the result of adding uniform, Poisson, and exponential noise to (a) exponentially decaying 
Wnt input, (b) linearly decaying input, and (c) exponential followed by linear decaying input. We add noise that is 
uniform on (0, 0.1), Poisson with mean 0.1, or exponential with a mean of 0.1. 



4. CONCLUSIONS 



It is difficult to draw conclusions from model predictions of response in the absence of noise (here used to 
model fluctuations in tumor response). While it is difficult to extrapolate models of tumor growth to in vivo 
outcomes due to inability to detect small tumors/^ the dramatic response of the ODE model to the addition 
of noise highlights the need to build better, stochastic models of tumor growth. The Markov model also shows 
how the response of the tumor fluctuates with time, which impacts on the predicted level of growth while 
undergoing treatment. Both the ODE and Markov models suggest a large tumor response to treatment is 
needed to dramatically improve the effectiveness of treatment. 

In modeling inter-cell signaling, omitting the noise terms results in lower than predicted /3-catenin levels. 
Given this, it is important to consider including noise when modeling in vivo /3-catenin; even small changes in 
Wnt-mediated /3-catenin expression result in clinically observable changes in tumor response.^' 

The models used are quite limited in terms of accurate modeling of the tumor growth and response, and the 
Wnt signaling. They make several broad assumptions, in particular the ODE model of tumor response assumes 
a homogeneous spherical mass and neglects the inner mass of necrotic cells surrounded by an irregular-shaped 
mass of rapidly dividing cells as dealt with in other models. Nonetheless they make useful predictions about 
cell-cell interactions and cell interactions with chemotherapy drugs. Quantitatively, adding noise to the tumor 
response in both the Markov and ODE models reduces the effectiveness in treatment; this is markedly shown 
in the ODE model. Adding noise to cell-cell interactions changes the levels of protein expression in the cell, 
which can then lead to an effect on tumor response. Building more accurate models would allow us to make 
more useful predictions of the effects of chemotherapy drugs. We propose further tests of these models in in 
vitro and in vivo settings, allowing further refinement of the models. 
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